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ABSTRACT 

Wc study the statistical properties of faint X-ray sources detected in the Chandra Bulge 
Field. The unprecedented sensitivity of the Chandra observations allows us to probe 
the population of faint Galactic X-ray sources down to luminosities L2-iokeV ~ 10^° 
erg s~^ at the Galactic Center distance. We show that the luminosity function of these 
CBF sources agrees well with the luminosity function of sources in the Solar vicinity 
(Sazonov et al. 2006). The cumulative luminosity density of sources detected in the 
CBF in the luminosity range 10'^°— 10^^ erg per unit stellar mass is i2-iokev/-^* = 
(1.7 ± 0.3) X 10^'' erg s~^ -^o^- Taking into account sources in the luminosity range 
10'^^ — 10'^'' erg s^^ from Sazonov et al. (2006), the cumulative luminosity density in the 
broad luminosity range 10^° - 10^'* erg s^^ becomes i2-iokcv/M* = (2.4±G.4) x 10^'^ 
erg per Mq. The majority of sources with the faintest luminosities should be 
active binary stars with hot coronae based on the available luminosity function of 
X-ray sources in the Solar environment. 



1 INTRODUCTION 

The hot interstellar plasma in galaxies (see e.g. Forman, 

Jones, & Tucker 1985) provides a valuable tool to probe 
different galactic phenomena, including energy and mass 
ejection from supernova explosions (e.g. Loewenstein & 
Mathews 1991; Mai'z-Apellaniz 2001; Strickland & Heckman 
2007), AGN activity (e.g. Padovani & Matteucci 1993; For- 
man et al. 2005), and the contribution of non-thcrmal pres- 
sure in hot atmosphores(Churazov ot al. 2010). However, 
before one can study the genuine X-ray emission of the hot 
interstellar plasma in early type galaxies, one must excise 
or account in some manner for the contribution of other 
constituents of galaxy emission. Revnivtsev et al. (2008) 
summarized the contributions from low mass X-ray bina- 
ries (see also Gilfanov 2004) and the fainter emission from 
the population of unresolved cataclysmic variables (CVs) 
and coronally active binary stars (ABs) (see Revnivtsev for 
a discussion of the origin of the faint diffuse emission) . 

The appearance of the latest generation of X-ray obser- 
vatories with fine angular resolution (e.g., as high as 0.6" 
from the Chandra X-ray Observatory) provides the possibil- 
ity to subtract/mask the contribution of the brightest X-ray 
binaries in nearby galax;ies. The contribution of fainter X- 
ray binaries is often estimated via extrapolation of the ob- 
served luminosity function of bright sources towards lower 
luminosities (e.g. Sarazin, Irwin, & Bregman 2001), which 
Gilfanov (2004) has shown to be a good estimator for this 
component. 

In the past, is has also been assumed that all emis- 



sion lines at energies <1-1.5 keV, visible in the spectra of 

galaxies, originate in a hot interstellar plasma (e.g. Sarazin, 
Irwin, & Bregman 2001; Irwin, Sarazin, & Bregman 2002; 
David et al. 2006). However, it is now clear (see (Revnivt- 
sev et al. 2006; Sazonov et al. 2006; Revnivtsev et al. 2007, 
2008, 2009, 2010)) that faint discrete galactic X-ray sources, 
namely cataclysmic variables and coronally active (usually 
binary) stars, provide an important contribution to the X- 
ray emission of galaxies (after subtraction of the contribu- 
tion from the brightest sources) and this should be taken into 
account. Moreover, as the X-ray emission from these stars 
naturally includes a contribution from optically thin plasma 
of moderate or low temperatures, their spectra contain ex- 
actly the same lines that are usually thought to originate in 
the interstellar plasma. 

Until now we have had very limited knowledge of the 
broad band luminosity function of the population of faint 
Galactic X-ray sources (however, there are studies of differ- 
ent sub-populations, see e.g. Pallavicini ct al. 1981; Worrall 
& Marshall 1983; Hiinsch et al. 1999; Schwope et al. 2002). 
Essentially, the best available broad band (10^^ erg < 
i2-iokev < 10^"* erg s~^) luminosity function was derived 
from all sky surveys of RXTE and ROSAT observatories 
by Sazonov et al. (2006). Deep Chandra observations of the 
peculiar Galactic Center region (Revnivtsev, Vikhlinin, & 
Sazonov 2007) added strong new support to the universality 
of the luminosity function over the entire Galaxy. 

Now, the new 1 Ms Chandra observation of the galactic 
bulge region - the Chandra bulge field (CBF) (see Revnivt- 
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sev et al. 2009, 2010) - allows us to probe the luminos- 
ity function of Galactic X-ray sources (mostly the old stel- 
lar population because sources are located in the Galactic 
bulge region, see Revnivtsev et al. 2010) to remarkably low 
luminosities. In particular, in the CBF, wc can probo to 
Lx ~ 10^° erg at the Galactic Center distance and thus 
to compare the derived luminosity function to that derived 
from the Solar vicinity (Sazonov et al. 2006) 

In this paper we explore the statistical properties of 
the discrete sources detected in the CBF and present their 
luminosity function and cumulative emissivity. 



2 DATA PREPARATION 
2.1 Chandra data 

Wc used all Chandra ACIS-I observations of the 
CBF (observation ID #5934, 6365, 6362, 9500, 9501, 
9502, 9503, 9504, 9505, 9854, 9855, 9892, 9893) with 
total effective exposure time ~ 898 ks. The Chan- 
dra data were reduced following the standard procedure 
fully described in Vikhlinin et al. (2005). The detec- 
tor background was modeled using the stowed dataset 
(http://cxc.harvard.edu/contrib/maxim/stowed). In this 
work we restrict our study to the circular region of 2.56 
arcmin-radius around l^^ = 0.113°, b'^ — —1.424°, close 
to the telescope optical axis (the region dubbed HRES in 
Revnivtsev et al. 2009), optimizing the angular resolution 
of the instrument and thus its sensitivity. 

To further enhance the sensitivity of the observation, we 
have searched for sources in the summed image of the HRES 
region in the broad 0.5-7 keV energy band. The sensitivity 
limit is /lim ~ lO"'^^ erg s"'^ cm~^. It is calculated for the 
energy spectrum of photons collected within the CBF in the 
energy band 0.5 7 kcV, approximated by a combination of 
models mekal with kT — 0.6 keV (abundance of elements 
was fixed at 0.1 of the Solar value) and bremss with kT = 7 
keV from the XSPEC package - the simplest model which 
allows us to reproduce the broad band properties of the en- 
ergy spectrum. Our model also includes interstellar absorp- 
tion with column density Nh ~ 7 x 10^^ cm~^. This value 
of interstellar absorption towards the CBF is adopted based 
on the results of Hong et al. (2009), measurements of Ay 
in Revnivtsev, Burenin, & Sazonov 2009; Revnivtsev et al. 
2010, and the Nh/Av ratio from Prcdchl & Schmitt 1995. 

Energy spectrum, collected from HRES is presented in 
Fig.l. Contribution of mekal spectral component to total 
0.5 7 kcV flux is 4.83 x 10"^^ erg s"^ cm~^ , contribution 
of bremss spectral component is 3.98 x lO^'^'^ erg s^'^ cm~^ 

The adopted sensitivity limit corresponds to a minimum 
detectable unabsorbed luminosity Lo.5-7 kev ^ 10^° erg s~^ 
at a source distance of 8 kpc, at which most of the Galac- 
tic objects in this field are expected to reside. In total, we 
have detected 473 sources with statistical significance > 4a 
(minimum number of counts per source is about 10). Statis- 
tical significance of detection was calculated with the help 
of rmmcrical simulations of the image taking into account 
the background and the exposure map. 

Due to the significant interstellar absorption of X-rays 
in the CBF and the great diversity of spectral shapes of 
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Figure 1. Energy specrum, collected by Chandra in HRES 
area. For conversion of Chandra counts into fluxes of 
sources we have adopted the simplest two component model 
(wabs* (mekal+bremss) model in XSPEC package, abundance of 

clciiicuts in mekal model is fixed at 0.1 of the Solar value), which 
allow us to reproduce the broad band properties of the energy 
spectrum. The spectral model was convolved with 100 eV Gaus- 
sian function to mimic a finite energy resolution of the instrument 

detected sources (e.g., the 2-7 keV/0.5-2 keV hardness ratio 

of cataclysmic variables and coronally active stars may differ 
by more than an order or magnitude), it is hard to make 
a definitive conversion of source count rates observed by 
Chandra into intrinsic fluxes and luminosities in the broad 
0.5-10 keV energy band. Therefore, despite the detection 
of sources in the broad energy band 0.5-7 keV, we present 
here only statistics of sources for the energy band 2-10 keV. 
Technically, the source fluxes were measured by Chandra 
in the energy band 2-7 keV and subsequently recalculated 
into the 2-10 keV energy band assuming the simple analytic 
form of the spectral fit described above, appropriate for the 
summed spectrum of all sources. We find that 1 count in 
the range 2-7 keV per image (898 ks efi'ectivc live time) 
corresponds to a source flux 4.6 x lO"^'^ erg s~^ cm~'^ in the 
energy band 2-10 keV. 

2.2 Stellar mass profile along the line of sight 

According to statistics of X-ray sources in the Solar vicinity 
(Sazonov et al. 2006), the majority of sources with lumi- 
nosities L2-iokev > 10^° erg s~^ are old binaries, including 
magnetic and non-magnetic cataclysmic variables, RS CVn 
systems, Algol-type binaries, W UMa-typa binaries. This 
allows us to safely assume that their density distribution is 
directly proportional to the stellar mass density distribution 
in the CBF. Therefore, it is useful to obtain the luminosity 
function of sources in the CBF normalized to the total mass 
of stars within the studied volume. 

The stellar density distribution along the line of sight in 
the CBF was previously studied in Revnivtsev et al. (2010). 



Luminosity function of faint Galactic sources in the CBF 3 



The main substructures of the Galaxy visible in the CBF are 
the Galactic disk, the Galactic bulge and the nuclear stellar 
disk (NSD). The majority of sources reside in the Galactic 
bulge and are located at ~8 kpc. In Revnivtsev et al. (2010), 
it was shown that the rms-width of the apparent magnitude 
distribution of red clump giants (RCG; considered to be al- 
most standard candles) is approximately 0.28'" (this rms 
scatter includes the intrinsic scatter of the RCG population, 
therefore the scatter caused by the distance variations alone 
is even smaller than 0.28'"). This means that the error in 
source luminosities, associated with the assumption that all 
sources in the CBF are located at 8 kpc, is less than ~ 30%. 

To determine the influence of these distance variations 
on the resulting luminosity function, and on our estimates 
of the luminosity density of sources, we have simulated a set 
of sources assuming that their intrinsic luminosity function 
has the shape dN/dL oc L^^ (which is close to what is actu- 
ally observed) and distributed the sources along the line of 
sight according to the stellar density distribution, described 
in Revnivtsev et al. (2010). The resulting luminosity func- 
tion does not change by more than 10% and the luminosity 
density changes by only ~ 5% in comparison to the case of 
constant source distances. 

According to the simple model of the stellar mass distri- 
bution (including the bulge, disk and the nuclear disk com- 
ponents) adopted in Revnivtsev et al. (2010), the total stel- 
lar mass angular surface density in CBF is ~ 4.6 x 10* Mq 
arcmin"^. We corrected for the fact that this stellar surface 
density corresponds to a ~ 13% higher NIR surface bright- 
ness than that measured by Spitzer in this area (Revnivtsev 
et al. 2010). Therefore in our subsequent analysis we adopt 
a stellar surface density 4.1 x 10* A/q arcmin"^ in the CBF 
and assume that the vast majority of RCGs lie at a distance 
8 kpc. This surface mass density is obtained via integra- 
tion of the adopted density profile up to further edge of the 
Galaxy, i.e. up to ~20 kpc distance from the Sun. Approxi- 
mately 85% of the all mass contained in small solid angles in 
the direction of CBF is located within 1 kpc of the Galactic 
Center. 



3 LUMINOSITY FUNCTION 

In calculating the number-flux function (or luminosity func- 
tion) of sources, we confined our analysis to the flux range 
where we expect no strong influence from incompleteness 
effects, i.e., where we detect more than ~ 10 counts from 
sources in the energy band used for source detection (0.5-7 
keV) and more than 5 counts from sources in the energy 
band 2-10 keV. In principle, it is possible to continue the 
luminosity function toward lower luminosities, but account- 
ing for the variety of spectral shapes at these luminosities 
(cf. soft spectra of coronally active stars and hard spectra of 
dwarf novae) introduces significant uncertainties - the same 
number of counts in total energy band, used for detection 
of sources, (0.5-7 keV) might correspond to significantly dif- 
ferent energy fiux in 2-10 keV energy band. 

The number-fiux function in its differential form, nor- 
malized to the mass of stars in the CBF area (4.1 X 10* Mq 
arcmin"^), is presented in Fig. 2. 

To account for the contribution of extragalactic sources 
(mostly AGNs) to the number-flux function detected in the 
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Figure 2. The 2—10 keV luminosity function of sources de- 
tected in the Chandra Bulge Field, corrected for the contribution 
of background active galactic nuclei. Numbers of sources, con- 
tributing to the presented bins of the luminosity function are 
1,7,15,31,76,103 from the highest to the lowest presented bin. 
For comparison we also present the luminosity function of faint 
sources obtained from the RXTE and ROSAT all sky surveys 
from Sazonov et al. (2006). The blue histogram shows the esti- 
mated (and subtracted) contribution of extragalactic sources in 
the CBF. The area unavailable to our survey due to the small 
number of sources per field of view is denoted by a dashed line 
dN/dlog L ~ 1. The luminosity of the brightest sources in our 
survey is marked by an arrow. 



CBF, we have adopted the results of Moretti et al. (2003). 
We calculated the anticipated number of AGNs in every lu- 
minosity bin, assuming the Moretti et al. analytic form of 
the AGN number flux function (in the energy band 2-10 
keV) . The AGN contribution has its own uncertainty mainly 
caused by either Poisson noise of the number of sources in 
the held, or by cosmic variance. According to simple es- 
timates of these uncertainties from the COSMOS survey 
(Cappelluti et al. 2007), in the small HRES field, the contri- 
bution from Poisson noise of the number of sources should 
dominate. We see that this noise is comparable with the 
Poisson noise of the number of Galactic sources in the field. 

The shape of the luminosity function is not clear, how- 
ever it can be approximated more or less adequately by a 
simple power law in luminosity interval 30.5 < logLx < 32: 
LdN/dlogL ^ A{L^ /lO'^^ erg/ sec)-"', where A « 1.1 x lO^'^ 
erg/s/Af0, and a = 0.2 ± 0.1 

The cumulative emissivity of Galactic sources as a func- 
tion of the limiting X-ray luminosity is shown in Table 1 and 
Fig. 3 (after correcting for AGNs as described above). Due to 
the relatively small survey volume, we cannot probe Galactic 
sources with luminosities higher than 10^^ erg s"'^ (increas- 
ing the solid angle of our study will not help here because 
the number of extragalactic sources will become a significant 
contributor to the source counts). Therefore, our estimates 
of the cumulative emissivity of sources lack the contribu- 
tion from higher luminosity sources (in particular, luminous 
cataclysmic variables). To compare our results with the lu- 
minosity density of sources determined in the Solar vicinity 
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Table 1. Cumulative luminosity densities of Galactic sources as 
a function of luminosity in the energy band 2-10 keV. Note that 
due to the limited survey volume we have only sources with lu- 
minosities less than 10^^ erg s~^. 
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Figure 3. Cumulative unit stellar mass emissivity in the 2-10 
keV energy band of the population of Galactic X-ray sources 
detected in the CBF (red area). The estimated contribution of 
AGNs was subtracted (see text for details). Since we do not probe 
sources with L2-l0keV > 10'^^ erg s~^ (indicated by the vertical 
dotted line), we have added the cumulative emissivity of sources 
at Z/2-lOkeV > 10^^ erg from Sazonov et al. (2006). The re- 
sulting shaded area with dashed perimeter can be compared with 
the corresponding result from Sazonov et al. 2006 (grey crosses). 



(Sazonov et al. 2006), we added the value of the cumula- 
tive emissivity measured for sources with lO^'^ erg s~^ < 
i2-iokeV < 10^* erg s~^ (shaded area with dashed perime- 
ter in Fig. 3). The resulting cumulative luminosity density 
of sources with 10^" erg s~^ < L2-iokcV < 10^'* erg s~^ is 
L2-iokov/M* = (2.4 ±0.4) X 10^7 



-lokcv < 10 erg s 
erg s """^per Mq. 



3.1 Contributions from various stellar populations 

It is practically impossible to determine the nature of all 
Chandra sources due to the extremely high spatial density of 
stars in the direction of the CBF (see e.g. images in Revnivt- 
sev et al. 2010). However, we can estimate the composition 
judging from our knowledge of the luminosity function of 
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Figure 4. Spectrum of the faintest sources (with 5-10 net counts 
per source over the whole exposure time, i.e. Lx ~ 10^" erg s~^), 
detected in the CBF (open circles). For comparison we present the 
spectrum of a typical bright active binary star II Peg (ASCA Ob- 
servatory data, observations Dec. 18-19, 1994), scaled to match 
the Chandra spectrum at energies 3-7 keV. The dashed curve 
shows a simple spectral model of a broken power law, which was 
fit to the spectrum and then absorbed with N[j = 7 X 10^^ cm~^, 
appropriate to the CBF. The curve has arbitrary normalization. 
It can be seen that the spectral shape of the faintest Chandra 
sources is compatible with that of bright ABs. 



different types of sources in the Solar vicinity. According to 
the space densities of different types of sources presented in 
Sazonov et al. (2006), we can anticipate approximately equal 
contributions from cataclysmic variables and from coronally 
active binaries to the cumulative 2-10 keV emisison in the 
luminosity range 10^^ erg s-i< L2 -lokev < 10 erg s , 
and a dominant contribution of ABs at energies below 2 keV. 

Unfortunately, measurements of the densities of cata- 
clysmic variables at lower luminosities (L2-iokeV < 10^^ 
erg s"'^) are more uncertain. Judging from estimates pre- 
sented in Pretorius et al. (2007) for CVs in the North Eclip- 
tic Pole region (Table 2 of this paper), their spatial density 
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cal stellar density of ~ 0.04 Mq pc~^ (Robin et al. 2003). 
The luminosity density of CVs can then be estimated to be 
~ 1.2 X lO^'' erg s"^ per Mq. This is just ~ 15% of the lu- 
minosity density of sources with 10'^" erg s~^ < L2-iokGV < 
10^^ erg s~^ in the CBF (~ (8 - 9) x lO^*' erg s'^ per 
Mq). This suggests that most of the detected sources in 
this luminosity range are active binary stars. 

Active binary stars produce softer X-ray emis- 
sion than do cataclysmic variables in general, so that 
io.5-2kov/i2-iokov > 1 for them (see e.g. Heinke et al. 
2005; Sazonov et al. 2006, see Giidel 2004 for review of X- 
ray emission of AB sources). This is of high importance for 
soft X-ray studies of other galaxies. In Fig. 4, we show the 
cumulative spectrum of the faintest sources detected in the 
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CBF (which is clearly softer than the summed spectrum of 
all detected sources) and compare it to the energy spectrum 
of the typical bright active binary star II Peg (see more de- 
tails on spectra of II Peg in Mewe et al. 1997). The dashed 
curve in the figure shows a simple broken power law fit to 
the spectrum of II Peg, arbitrarily scaled and absorbed with 
interstellar column density Nh = 7 x 10^^ cm~^, appropri- 
ate for the CBF. We see that the spectrum of the faintest 
CBF sources is compatible with that of AB stars of the ap- 
propriate luminosity. However, the considerable absorption 
at energies < 2 keV prevents seeing the full luminosity den- 
sity of such soft sources. Without interstellar absorption, 
the CBF brightness at energies 0.5 — 2 keV would be several 
times higher. 



4 SUMMARY 

We have analysed statistical properties of X-ray sources de- 
tected in deep observations of the Chandra bulge field. The 
results of this work can be summarized as follows: 

• The luminosity function of detected X-ray sources in 
the luminosity range probed by our observations (30 < 
logL2-iokeV < 32) and corrected for the contribution of 
extragalactic sources (AGNs), is relatively flat {L^dN/dL ~ 
constant) and consistent with that determined from a sam- 
ple of sources in the Solar vicinity by Sazonov et al. (2006). 

• The cumulative luminosity density of sources in the 
luminosity range 30 < logL2-i()koV < 32 is (1.7 ± 
0.3) X lO'^^ erg s^^per Mq. Accounting for the contribu- 
tion from sources with logL2_iokGV > 32 from Sazonov et 
al. (2006), the cumulative luminosity density of sources with 
logL2-iokev > 30 is (2.4 ± 0.4) x lO^'' erg s"^ per Mq. 

• The majority of the faintest sources should be active bi- 
naries with hot coronae, which makes their cumulative emis- 
sion virtually indistinguishable from that of the hot inter- 
stellar plasma typically associated with old stellar popula- 
tions. This indicates the importance of carefully accounting 
for the emission of such sources, based on detailed knowl- 
edge of their statistical properties and not only on the shape 
of the energy spectrum. 
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